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ABSTRACT 

We present the analysis of the color-magnitude relation (CMR) for a sample of 57 X-ray detected 
Abell clusters within the redshift interval 0.02 < z < 0.18. We use the B — R vs R color-magnitude 
plane to establish that the CMR is present in all our low-redshift clusters and can be parameterized 
by a single straight line. We find that the CMRs for this large cluster sample of different richness and 
cluster types are consistent with having universal properties. The k-corrected color of the individual 
CMRs in the sample at a fixed absolute magnitude have a small intrinsic dispersion of ^ 0.05 mag. 
The slope of the CMR is consistent with being the same for all clusters, with the variations entirely 
accountable by filter band shifting effects. We determine the mean of the dispersion of the 57 CMRs 
to be 0.074 mag, with a small rms scatter of 0.026 mag. However, a modest amount of the dispersion 
arises from photometric measurement errors and possible background cluster superpositions; and the 
derived mean dispersion is an upper limit. Models which explain the CMR in terms of metallicity 
and passive evolution can naturally reproduce the observed behavior of the CMR in this paper. The 
observed properties of the CMR are consistent with models in which the last episode of significant 
star formation in cluster early-type galaxies occurred significantly more than ~ 3 Gyr ago, and that 
the core set of early-type galaxies in clusters were formed more than 7 Gyr ago. 

The universality of the CMR provides us with an important tool for cluster detection and redshift 
estimation. A very accurate photometric cluster redshift estimator can be devised based on the 
apparent color-shift of the CMR due to redshift. This calibrator has the additional advantage of 
being very efficient since only two bands are needed. An empirically calibrated redshift estimator 
based on the B — R color of the CMR for clusters with z < 0.20 produces an accuracy of Az ~ 0.010. 
Background clusters, typically at z ~ 0.25 and previously unknown, are found in this survey in the 
color-magnitude diagrams as secondary CMRs to the red of the target cluster CMRs. We also find clear 
cases of apparent X-ray substructure which are due to these cluster superpositions. This suggests that 
X-ray observations of clusters are also subject to a significant amount of projection contaminations. 
Subject headings: Galaxies: clusters: general — galaxies: elliptical and lenticular, cDs — galaxy: 
formation — galaxies: photometry 



1. INTRODUCTION 

Early-type galaxies are the domin ant population 
in low-redshift galaxy cluste rs (e.g., iDressled 119841: 
iKormendv fc Diorgovskil [19891) . When compared with 
other galaxy types, early-type galaxies can be labeled 
as a homogeneous class in terms of the dynamical 
structure, stellar populations, and the gas and dust 
content. The remarkable uniformity of the properties 
of early-type galaxies in clusters and the tightness 
of their fundamental plane and its projections sug- 
gest that these galaxies are coeval, had their last 
episode of strong star formation at high redshift 
(z > 2), and have evolved passively since the n (e.g., 
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As larger samples of early-type galaxies become avail- 
able, the idea that giant spheroidal galaxies had their 
major starburst at high redshifts, regardless of the 
envir onment, becomes more popular (e.g., iPeeblesI 

Hoi. 

One of the properties of early-type galaxies that can 
be studied most straightforwardly is the color-ma gnitude 
relation (CMR). This effect was discovered bv IBaumI 
ljl959|) who noted that field elliptical colors become red- 
der as the galaxies become brighter. 'Baum's study also 
revealed that the average colors of ellipticals are much 
redder than those of glob ular clusters, and concluded, 
contrary to iBaadj l)1944j) . that elhptical galaxies were 
dominated by old Population I stars. To arrive at this 
conclusion, Baum generated a simple population syn- 
thesis model that was based on the, then novel, idea 
of progressive metal enrichm ent of a population by suc- 
cessive gener ations of stars ijFowler fc GreensteirJll95^ 
lStruvelll956D . Hence, elliptical galaxies may have ex- 
perienced an extended period of very efficient star for- 
mation during the early stages of their formation. This 
idea is now accepted, but still debatable in terms of our 
inability to disentangle the effects of age and metallic- 
ity in present day ellipticals ( e.g., iWorthev et al.llT995t 
iFerreras. Cha,r]ot. fc Si]kllT99l. 

A simple straight line fit can describe the CMR 
for elliptical galaxies in an interval of about 8 mag- 
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nitude s in Virgo (e.g., iSandaed 19721) and Coma 
(e.M.. iThomnson CregorylTmM |L6pez-Crnz| dH 
iSecker. Harris, fc Plummerf Il998|) . This effect IS re- 
markable because it covers a range in luminosity of 
about 10^, suggesting that within this range galax- 
ies have shared a similar evolutionary process. How- 
ever, this is not the only way to get such an effect. 
For instance, the CMR can be reproduced by demand- 
ing younger galaxies to be m ore metal rich than older 
galax ies of similar luminosity ('Ferr eras. Chariot, fc Sillj 
Il999|) . The potential of the CMR as a tool f o r cos- 
mological studies was recognized by iSandagel l|1972f) 
and was characterized for different galax y environments 
(^isvanathan & Sandag3 ll977t [Sandagc fc VisvanathanI 
Il978albf) . The CMR's application as a distance indicator 
has been suggested by Visvanatha. n fc San dagc ( 1977 .) • 
The dispersion in the CMR, which is very small, has been 
quantified for Virgo and Coma bv lBower! Lucev. fc ElTi^ 

(^992, hereafter BLE) and for high-redshift clus- 

ters (e.g., Elhs et al. 1997; Van Dokkum et al.' '199^; 
[St anford TEiscn hardt , & Dickinson 1998; Glad ders et aL , 
Il998: Blakcslc e et al.ll2(X)3)) . These studies have used the 
dispersion to constrain the epoch of galaxy formation and 
the way they have evolved. 

The CMR was modeled as a metallicity ef- 

fect bv [Arimoto fc Yoshii l|1987t) . The mo dels of 
IArim9t o fc Yoshiil incorporated the ideas of iLarsonl 
l(T974afa) in whkh early-type galaxies are formed by 
dissipative monolithic collapse and the star forma- 
tion regulated by supernova driven winds. As in 
[Bower et al. (1998J, we will call Larson's mode l the 
classical model. More re cently, iKodamal (^1997^ and 
iKodama fc Arim g^ l^^I l- hereaft er KA97) have revised 
the models of f^imotofc Yoshiil (^87) and have in- 
cluded improved libraries for stellar evolution, incorpo- 
rating the effects of the variation of metallicity with 
galaxy mass. Alternatively, t he sem i-analytical model- 
ing of iKauffmann fc Chariot! l|1998fl . which incorporate 
the hierarchical galaxy formation scenario, can also re- 
produce the CMR for early-type galaxies by imposing 
a very early epoch for the major mergers. Much effort 
has been devoted to improving the models and the ob- 
servations; however, it is still difficult to decide whether 
early-type galaxies for med according to t he hierarchical 
or the classical model (jBower et al.lll99^ and references 
therein). 

If cluster early-type galaxies were formed at high red- 
shift, are coeval and evolving passively, then it should be 
expected that the CMR is universal at the current epoch. 
Universality should be defined in the present context as 
the validity for all clusters of the same rest-frame pa- 
rameters: slope, intercept and dispersio n. The works of 
Sandage fc VisvanathanllVisvanathan fc~andagel 119771 
iSandage fc Visvanathanlll978alhD and BLE made a case 
for the universality of the CMR for early-type galax- 
ies, albeit based on a few local clusters. There are 
other studies that appear to agree with these results 
/e.g. IButcher fc OemleHrmsl inixon. Codwiu. fc Pea,chl 



[T989t IMetcalfe et al][l99l ICarilH et al.|[l99fil and ref- 
erences therein). However, some apparent cluster- 
to-cluster variations of the CMR have also been re- 
port ed, but the evidence is not c o mp elling. For exam- 
ple, lAaronson. Persson. fc Frogell l)1981il reported that 
the U — V and V — K colors for galaxies in Coma 



were bluer than those for galaxies in Virgo; this re- 
sult was not confirmed by BLE. Drastic changes in 
the CMR's slope for gala xies within the same cluster 
have been reported (,Metc alfe et al.l[T9 '94. and references 
th erein). The model of gala ctic cannibalism proposed 
bv lHausman fc Ostriked l|1978,l hinted at an explanation 
for this phenomenon. However, most of the studies at 
differ ent redshifts, even as high as z ~ 0.5 llEllis et al 
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1998 : lGladders et al .' 1998), have shown little indications 



of such an effect. The studies mentioned above utilized 
relatively small cluster samples and accurate photometry 
was limited mainly to the few nearest clusters. In this 
paper we show that the CMR in a 57-cluster sample at 
low redshifts, with different richness, cluster types and 
X-ray luminosities, are consistent with being universal. 

This paper is the second of a series resulting from 
a large multicolor imaging survey of low-redshift Abell 
clusters. The paper is organized as follows. In §2 we 
present the sample selection criteria. In §3 we briefly de- 
scribe the observations and photometric reduction pro- 
cedures. Section 4 describes the CMR parameterization 
and the error estimate techniques; while §5 examines the 
dispersion of the CMR. We present a discussion of the 
results, in particular regarding the universality and the 
evolutionary models, in §6. In §7 we introduce some 
observational applications of the CMR based on the re- 
sults of this study, including using the CMR as a redshift 
indicator and a cluster finding tool. Finally in §8 we 
summarize the conclusions. More details regarding sam- 
ple selection, observations, image pr eprocessing, catalog s 
and finding chart s can be found in iLopez-Crud l)1997ji . 
iBarkhousg ({200.'^ . and Barkhouse, Lopez-Cruz, fc Yee 
(2004a, in preparation. Paper I of this series). A detailed 
discussion of the luminosity function of cluster galaxies 
using this survey can be found in Paper HI (Barkhouse, 
Yee, fc Lopez-Cruz 2004b, in preparation). Paper IV de- 
scribes the characteristics of various color-selected clus- 
ter galaxies (Barkhouse, Yee, fc Lopez-Cruz 2004c, in 
preparation). Recent cosmological observations indicate 
that the best world model agrees with a flat A-dominated 
universe that is characterized by r jm — 0.3, — 0. 7 
and Ho ~ 70 km s'^ Mpc'^ (e.g.. iSnergel et al.ll200l . 
However, considering that the effects of dark energy and 
curvature are negligible at low redshifts {z < 0.2) and in 
order to allow direct comparisons with previous studies; 
we, then, for convenience have set, unless otherwise indi- 
cated, Ho = 5O/150 km s~^ Mpc""'^ and qo = throughout 
this paper. 

2. THE SAMPLE AND OBSERVATIONS 

The sample contains Abell clusters selected mainly 
from a compilation o f bright X-ray cl usters from Ein- 
stein's IPC made bv iJones fc FormanI (^99). The se- 
lection of clusters identified by independent techniques 
(i.e., optical and X-ray) reduces the probability of con- 
fusing real, gravitationally bound clusters, with appar- 
ent galaxy over-densities due to projection effects. The 
initial sample was defined under the following selection 
criteria: 1) the clusters should be at high galactic lati- 
tude, |6| > 30°; 2) their redshifts should lie within the 
range 0.04 < z < 0.20; 3) the Abell richness class (ARC) 
should, preferably, be greater than 0; and 4) the declina- 
tion 5 > 20°. We attempt to follow these criteria strictly; 
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however, some ARC=0 clusters were included due to ob- 
servational constraints, such as the lack of suitable clus- 
ters at certain right ascensions during the observations. 
The ARC=0 clusters were not selected at random, but 
on the appearance of the X-ray emission and the Bautz- 
Morgan type that indica te a higher richness class. See 
lYee fc L6nez~C>^ itTOQfll) for a discussion of the richness 
of the cluster sample. 

This sample includes 47 clusters of galaxies observed 
in Kron-Cousins B, R, and / at KPNO with the 0.9m 
telescope using the 2048 x 2048 pixels T2KA CCD 
(IT.6nez-Cru7ll1997HYee T.onez-Cru^rm flflt IT.onez-Cnl^ 
I2001L hereafter the LOCOS sample). In addition, two 
clusters imaged in B and R were included from Brown 
(1997) using the same instrumental setup and satisfy- 
ing our selection criteria. We note that only the B 
and R data are considered in this paper. The field cov- 
ered by this combination of a small telescope plus a 2K 
CCD is 23.2' X 23.2' with a scale of 0.68"/pixel, i.e., 

l.bh^Q Mpc at z = 0.04 and ^ 6fe^n Mpc at z = 0-2. 

A sub-sample of eight clusters from lBarkhoiis3 l)2003f l 
are included to complement our 49 cluster sample by 
covering the low redshift interval from 0.02 < z < 0.04. 
The data were obtained at KPNO with the 0.9m tele- 
scope using the 8k Mosaic camera (8192 x 8192 pix- 
els). This telescope/detector combination yields a scale 
of 0.423"/pixel, with a field-of-view of one square degree 
(~ 2h^^ Mpc at z = 0.02 and ~ ih^^l Mpc at z = 0.04). 
The clusters from this sample were selected using the 
previous criteria except that ARC=0 clusters were not 
preferentially excluded, although all clusters in the sam- 
ple are detected in X-rays and have a prominent CMR. 

The integration times for our 57 cluster sample varied 
from 250 to 9900 seconds, depending on the filter and the 
redshift of the cluster. The pho t ometr ic calibrations were 
done using stars from lLandold 1)1992(1 . Control fields are 
also an integral part of this survey. We observed in both 
R and B filters a total of 6 control fields chosen at random 
positions on the sky at least 5° away from the clusters 
in the sample. These control fields were observed using 
the T2KA CCD and the Mosaic camera to a comparable 
depth and reduced in the same manner as the cluster 
data. All observations included in this study were carried 
out during 1992-93 and 1996-98. 

Table 1 lists the sample with the cluster redshifts. 
Coma (A1656) was included in the LOCOS sample de- 
spite its low redshift (z — 0.0232). The main reason for 
its inclusion is the vast amount of data for Coma avail- 
able in the literature, which allows us to make direct 
comparisons of our results with others. (We note that 
only R images of Coma were obtained with the Mosaic 
camera and hence are not considered here.) 

3. PHOTOMETRIC REDUCTIONS 

The preprocessing of the images was done using IRAF. 
The photometric reduction was carried out using the pro- 
gram PPP(Picture Processing Package. r^ll99H) . which 
includes algorithms for performing automatic object find- 
ing, star/galaxy classification and total magnitude de- 
termination. We also exploited a series of improvements 
to RPR described in fVce. Ellingson. & Carlbcre (1996) 
which decrease the detection of false objects and allow 
star/galaxy classification in images with variable point- 
spread function (PSF). We describe the reduction proce- 



dures briefly below; full details are given in Paper I. 

3.1. Photometry 

The object list for each cluster is compiled from the 
R frames. The R frames are chosen because they are 
deeper than the images from the other filters. On av- 
erage, about 3000 objects are detected in each LOCOS 
fields and about 25,000 objects in the mosaic fields. Var- 
ious steps are taken to minimize spurious object detec- 
tions due to bleeding columns, bright star halos, artificial 
satellite tracks and other cosmetic effects. 

PRR uses the curve of growth of counts to determine 
total magnitudes of identified objects. A sequence of 30 
concentric circular apertures is used, starting from a di- 
ameter of ^ 2" (three pixels), approximately the diam- 
eter of the seeing disk. The apertures have a maximum 
diameter of 20 - 24", in steps of 0.68" or 0.423" (i.e., one 
pixel). An optimal aperture size for each object is deter- 
mined based on the shape of th e growth curve, using the 
criteria described in lYed l|199H) . 

Determining the photometry for galaxies in clusters is 
more complicated than for field galaxies at comparable 
magnitude limits. One of these problems is the large 
range of galaxy sizes, which extends from a few kpcs 
for dwarf galaxies to ~ l^so^ Mpc for cD galaxies. We 
find that for bright galaxies whose apparent R magni- 
tude is brighter than about 18.5, an aperture of 20" is 
insufficient for the determination of the total magnitude. 
Hence, in a second iteration, the growth curves of all ob- 
jects brighter than R — 18.5 and classified as galaxies 
are extended to a maximum aperture of 81 — 120" with a 
coarser radial increment, allowing the total light of even 
the largest non-cD galaxies to be measured at z ^ 0.02 

The presence of close neighbors and crowding is a prob- 
lem that is more acute in cluster galaxy photometry. PRR 
deals with this by automatically masking all neighbors 
within twice the radius of the maximum allowable aper- 
ture from the object being considered (see Ycc 1991). An 
additional problem which may adversely affect the pho- 
tometric and color accuracy of cluster galaxies situated 
near the cluster core is the presence of a cD galaxy that, 
by virtue of its large size, may engulf many nearby galax- 
ies in its envelope. To solve this problem, the photome- 
try of galaxies near the core is carried out after the cD 
and bright early-type galaxies have been removed usin g 
profile modeling techniques developed bv IBrowi] 1)1997(1 . 

3.2. Color Determination 

The colors for the galaxies are determined using fixed 
apertures on the images of each filter, sampling identical 
regions of galaxies in different filters. Due to the large 
range of cluster galaxy sizes and the sampled redshifts 
range (0.02 < z < 0.18) we adopt the following scheme. 
For bright galaxies (i? < 17.5) at z < 0.06, a maximum 
physical aperture of O.OII/i^q^ Mpc is used. This aper- 
ture varies in angular size from 19" for z — 0.04, to 6" 
for z = 0.06. For bright objects at larger redshifts the 
maximum allowable aperture is fixed at 6", this corre- 
sponds to 0.025/i^Q^ Mpc at z — 0.18, the largest cluster 
redshift in the sample. This change in aperture for the 
objects at redshift larger than z = 0.06 is introduced 
because the angular size of an O.Oll/igQ^ Mpc aperture 
becomes too small at those redshifts. Using too small an 
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aperture introduces systematic errors due to seeing and 
variations of the PSF; hence, a minimal color aperture of 
~ 3x FWHM is used to avoid these effects (the average 
seeing measured in the cluster images is ~ 1.5"). Overall 
the internal accuracy in the color determinations should 
be about 0.005 magnitudes in _B — i? for bright objects. 
For faint objects, the physical size is often smaller than 
0.011/i^Q^ Mpc; hence, the smallest optimal aperture from 
the available filters as determined by the growth curve 
profile is used. The errors for faint objects can be as 
large as 0.5 magnitudes in B — R. 

The approach of using relatively small aperture sizes 
for color determination increases the accuracy, because 
only the central region of each object is used. This, how- 
ever, comes at the expense of assuming no color gradi- 
ents in the galaxies. The effects of color gradients in 
B ~ R are very small for early-type galaxies; color gradi- 
ents of only A(. B-i?)/Alogr = -0.q9±0.02 rnag per dex 
in radius fe.g.. IPeletier et al.llT990t IWu et al.ll2004D are 
present in elliptical galaxies. However, larger color gra- 
dients are found in late- typ e galaxies (e.g., de Jong 199_S 
iGadotti fc dos Aniosl200l() : galaxies could be bluer by as 
much as B — R = 0.8 mag between t he nuclear and the 
outer regions fsee Ta ylor et al.l 120031 for an alternative 
viewpoint). However, those color-gradients do not se- 
riously affect our color determinations, since these very 
late-type galaxies are not a significant population in rich 
clusters of galaxies for our redshift range. 

We note that the total magnitude of a galaxy is deter- 
mined using the growth curve from the R image, while 
the total magnitudes in the B and / images are deter- 
mined usin g the color differences with respect to the R 
images fsee T^ITMl for more details). 

3.3. Star/ Galaxy Classification 

Star-galaxy classification is a very important issue in 
wide-field galaxy photometry, because the foreground 
stellar contribution relative to galaxies is large at the 
bright end and has large variations from field to field. 
PPP uses a classifier that is based on the comparison of 
the growth curve of a given object to that of a reference 
PSF. The reference PSF is generated as the average of the 
growth curves of high signal-to-noise ratio, non-saturated 
stars within the frame. The classifier measures the "com- 
pactness" of the object by effectively comparing the ratio 
of the fluxes of inner and outer parts of an object with 
respect to the reference PSF. 

An additional problem with wide-field imaging is the 
PSF variation across the frame, which, although not se- 
vere, is clearly detec table in these data, especially for 
the mosaic images. lYee. Ellingson. fc Carlberel l|199fif) 
describe a procedure of using local PSFs for star/galaxy 
classification to compensate for this effect. Local PSFs 
generated from a set of bright stars, ranging in number 
from 30 to 200 per frame and distributed over the whole 
image, are used for classification in these images. The 
detailed procedure, including visual verifications and re- 
liability, is discussed in Paper I. 

3.4. Calibration to the Kron-Cousins Standard System 

Instrumental magnitudes are calibrated to t he Kron- 
Cousin s system by observing standard stars from lLandoltl 
1)1992(1 . Due to the large field, up to 45 standard stars 



can be accommodated in a single frame. The color prop- 
erties of the standard stars cover a large color range that 
encompasses those of elliptical and spiral galaxies. 

We observed standard star fields in each filter ^ three 
times throughout the night. The standard stars are mea- 
sured using a fixed aperture of 30 pixels for the LOCOS 
frames and 32 pixels for the mosaic data. These aperture 
sizes are selected as being the most stable after mea- 
suring the magnitudes using a series of diameters. We 
adopt the average extinction coefficients for KPNO, and 
fitted for the zero points and color terms. The airmass 
terms for the LOCOS data are held fixed to the values 
-0.270, -0.100 and -0.040 (in units of magnitudes per 
air mass) for B, R, and /, respectively. For the mosaic 
observations, the adopted average extinction coefficients 
are fixed at -0.243 for B and -0.097 for R. Nightly so- 
lutions for the remaining coefficients are obtained. The 
rms in the residuals of individual fittings is in the range 
0.020-0.040 mag, which is comparable to the night-to- 
night scatter in the zero points. This can be considered 
as the systematic calibration uncertainty of the data. 

3.5. Completeness and Final Catalogs 

The final galaxy catalogs are generated using all the 
information and corrections derived in the previous sec- 
tions. For data obtained under non-photometric condi- 
tions, single cluster images were obtained during photo- 
metric nights in order to calibrate the photometry (three 
clusters in total). A final step involves the determination 
of the completeness limit. A fiducial 5a detection limit 
is determined for each field by calculating the magnitude 
of a stellar object with brightness equivalent to having 
a S/N = 5 in an aperture of 2". This is done by scal- 
ing a bright star in the field to the 5a level. However, 
the 5a limit is fainter than the peak of the galaxy count 
curve and hence is below the 100% completeness limit 
for galaxies. A conservative 100% completeness limit is 
in general reached at 0.6 to 1.0 magnitudes brighter than 
the 5a detection. Seeft^ce (1991) for a detailed discus- 
sion of the completeness limit relative to the 5a detection 
limit. 

Each cluster's final catalog contains the object iden- 
tification number, pixel positions, the celestial coordi- 
nate offsets with respect to the brightest cluster galaxy 
(BCG), apparent magnitudes B, R and / (/ magnitudes 
are only available for 47 clusters), and their respective 
errors, object classification, and R — I and B — R colors. 
A more detailed presentation of the catalogs and gray 
scale images of the clusters is given in Paper I, and can 
be accessed in the near future from a Web site. 

3.6. Galactic Absorption Correction 

In order to compare the galaxy colors in a consistent 
manner, we have to correct for the extinction produced 
by our own galax y. The corrections for t he filters used 
in this study are l|Postman fc Laueilll995D : 

Ab = 4.05E{B - V), 
Ar^2.35E{B -V). 

The values of the g alactic extinction coefficie nts can be 
calculated from the iBurstein fc HeilesI l)1982D maps, us- 
ing the reported E{B — V) values, o r directly from the 
Ab tabulations for bright galaxies l|Burstein fc Heilel 
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119841) with coordinates in the vicinity of our pointed 
observations using NED. Alternatively, extinction maps 
based on dust infrared emission have been provided by 
[SchlcEcl. Finkbciner. & Davis (1998). We found that Ab 
values, for th e cluster in ou r sample, determined by 
ISchleeel . Finkbeiner. fc Da^^ j are systematically larger 
than those of [Burstcin & Heiled , but in most cases the 
difference is not larger than 0.1 mag. In order to make 
direct comparison with previous works, we have adopted 
iBurstein fc HeilesI determinations. We are aware that 
both schemes have their own limitations and biases: the 
difference and uncertainties between these two extinction 
schemes are, therefore, considered as a systematic error. 
In Paper I we provide the extinction values used for each 
cluster. 

3.7. The Color- Magnitude Diagram 

We plot the B — R color versus i?-band magnitude to 
create the color- magnitude diagrams (CMD). Both the 
magnitudes and colors have been corrected for galactic 
extinction. The entire set of cluster CMDs can be found 
in Paper I. As examples, Figure 1 shows the CMDs for 
six of the clusters, demonstrating the various properties 
of the CMR. F ollowi ng the terminology introduced in 
IMetcalfe et all lll99l we can readily distinguish three 
zones: the "blue" zone, which is populated by late-type 
galaxies both in the cluster and in the field, and early- 
type galaxies with redshifts lower than that of the clus- 
ter; the "sequence" , which contains galaxies that have 
a high probability of being e arly-type cluster members 
at the cluster reds hift (e . g.. IPixon. Godwin, fc PeachI 
19891 iBiviano et alJ 119951 iThompson fc Gregorvl 119931: 
Yee et al.l 119961 IVazdekis et al.ll200H) : and finallv. the 
"red" zone, which is populated by galaxies which have 
redshifts higher than that of the cluster. 

4. THE PARAMETRIZATION OF THE CMR 

Different approaches have been applied to quantify 
the CMR. In most cases, the analysis has only been 
applied to objects with confir med early-type morphol- 
ogy and membership (e.g., BLE: lVisvanathan fc Sandagd 
[1977). Because the presence of outliers is minimized by 
membership and morphology information, simple least- 
squares methods can be applied in these cases. In the 
present study we do not have redshift information for 
the majority of the cluster galaxies nor galaxy morpho- 
logical classifications. Nevertheless, as illustrated by the 
CMDs, the CMR stands out, and it can be easily traced 
visually. Robust methods can be applied to provide a 
statistical fit that can be resistant against the influence 
of outliers. A few schemes include those used by BLE, 

fetcalfc et al. (1994), Ellis et al. (1997), Gladdcrs^Tal] 
998), and fPimbblet et alJ l(2002i) . Below we describe 
the robust fitting technique used in this paper. 

4.1. Robust Regression Based on the Biweight Method 

The red sequence of early-type galaxies can be pa- 
rameterized using a simple straight line. Fitting this 
line to the red sequence of the CMR is not an unam- 
biguous procedure. In practice, the presence of a sig- 
nificant number of outliers, such as background galax- 
ies, tends to skew the derived solution. Simple least- 
squares methods are ineffective in this situation and more 



"robus t" estimators, suc h as the Least Ab solute Devi- 
ation press et alJ 11992^1 . have been used l|L6pez-Crud 
Il997|) . After numerous iterations using a variety of 
strategies, a technique utilizing the biweight method 
llBeers. Flvnn. fc Gebhardtl Il990) was adopted. This 
method is found to be very robust and not greatly in- 
fluenced by the presence of outliers. A variant of this 
method has been implem ented by sev eral gr oups to ef- 
fectively fit the CMR (e.g..lTerlevich." Caldwci rfc Bowed 
l20()ltlPimbblet et a l. '2002\! 

The procedure used in this study to successfully fit the 
CMR is summarized as follows: An initial cut is made to 
each cluster catalog to cull galaxies further than 0.2 7-200 
from the adopted cluster center, which is typically the cD 
or brightest cluster galaxy. The value of r2oo , the radius 
within which the average density is 200 times the critical 
density, is derived p hotometrically for each cluster based 
on cluster richness l|Yee fc Ellings onM2003^. see paper III 
for a detailed discussio n), and is a pproximately equal to 
the virial radius (e.g.. iCole fc Lacevlll99flll . The radius 
criterion provides a less biased measurement of the CMR 
since our cluster images cover a large range in linear size. 
By using the same dynamical radius for each cluster we 
minimize the influence of using different cluster-centric 
radii, which can potentially affect our CMR measure- 
ments (e.g., Pimbblet et al. 2002). The value of 0.2r2oo 
was chosen since it is the maximum radius in which all 57 
clusters are fully observed. For this initial cut of the clus- 
ter catalog, galaxies fainter than ^ 3 magnitudes brighter 
than the completeness limit are removed in order to help 
reduce the presence of outliers. Once the modified clus- 
ter catalog is generated, an estimate of the slope and 
intercept is made by visually inspecting the CMR. A set 
of deviations given by di = yi — (a + hx), where Xi and yi 
are the observed values of R and B — R for each galaxy, 
b is the slope, and a is the intercept, is calculated for 
each galaxy. A biweight function is then used to calcu- 
late the location and scale. In this context, the location 
and scale are statistical quantities that describe where 
the data are mostly found (e.g., median) and the spread 
in the data about this value (e.g., dispersion), for the set 
of deviations (Beers et al. 1990). This function is given 

where u={di- M)/{c- MAD), M is the median, MAD 
is the median absolute deviation from the median, and c 
is the "tuning" constant (taken to be 9.0, see Beers et al. 
1990). The value of the slope and intercept are than 
incremented by a small amount e and a new set of devia- 
tions constructed. The location and scale are determined 
for this new set, and the process is repeated for a range 
of slope and intercept values that bracket the estimated 
best-fit values. The grid of slope and intercept values are 
searched to determine the best-fit value of a and b that 
minimizes the location and scale of the deviations. 

4.2. Error Estimate and Outlier Rejection 

We apply the n on-parametric bootstrap method 
l|Efron fc Tibshiranl 1986: Babu fc Rao 1993) to derive 
the errors. iBabiifc Singhl l(l98.^ have analytically es- 
timated that N n log^ n (where N is the number of 
resamples to generate and n is the number of elements 
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in the original data set) resamples give a good approxi- 
mation of the underlying density distribution. We adopt 
the median values for the slope and intercept derived 
from N = n log^ n bootstrap resamplings as the solution 
for the fit. These values are slightly different from the 
ones derived from the direct fit of the actual data. The 
standard errors for the median slope and intercept are 
the rms of these quantities from the bootstrap samples. 

In some clusters the standard errors of the slope are 
large, with the uncertainty comparable to the slope it- 
self. We identify three reasons for this occurrence: a) the 
presence of a second CMR produced by a serendipitous 
cluster or group of galaxies at a higher redshift; b) when 
the cluster is poor and there are gaps (in R magnitude) 
in the CMR; and c) for the highcr-rcdshift clusters where 
the clusters show indications of the Butcher-Oemler ef- 
fect l|Butcher fc OemledlTop) . 

For those clusters with large error estimates, a second 
pass with outlier rejection is performed which allows the 
fitting to move closer to the line defined by the CMR and 
produce more reasonable error estimates. Points that are 
2(7 away from the main fit are rejected and the procedure 
is repeated. Three cases can be recognized in our second 
pass: a) no effect at all, b) producing a similar fit as the 
first pass but with a smaller error estimate, or c) passing 
from one unique fit to another one. By inspection, we 
normally recognize that the initial fit is already a good 
fit to the CMR. The second pass is introduced to refine 
the solution and reduce the estimated standard errors. A 
three sigma rejection criterion usually results in case (a). 
Rejection at the one sigma level invariably ends in case 
(c). The optimal approach is obtained using a two sigma 
rejection procedure that results in a new fit consistent 
with the previous iteration, but with smaller computed 
errors, i.e., case (b). In some cases where the CMR is 
poorly populated, the errors become larger after applying 
the rejection criteria, a consequence of a large change in 
the fit. We note that only three clusters required a second 
pass with outlier rejection. The method outlined above 
was found to produce an accurate characterization of the 
CMR, as can be seen by visual inspection of the color- 
magnitude diagrams in Figure 1. This procedure was 
also used to success fully fit the CMR of A2271 and A2657 
from the sample of iLopez-Crud lll997j) . which could not 
be achieved by utilizing the Least Absolute Deviation 
jPress et al...l99 2) method used in that study. 

4.3. Results 

The resultant parameters of the CMR fits are given in 
columns 3 and 4 of Table 1, which tabulate the slope and 
intercept and their uncertainties for each cluster. The fit 
of the CMR is shown in the example CMDs in Figure 1 
as a dashed line with the fit equation shown. The plot 
of the fits of all the clusters can be found in Paper I. 

5. THE DISPERSION OF THE CMR 

The dispersion of the CMR is a useful diagnostic for 
the star formation history of early-type galaxies in clus- 
ters (e.g., BLE). In order to provide a robust measure 
of the dispersion of the CMR, we create for each cluster 
(using an identical radius cut as that used in determining 
the CMR fit parameters, §4.1) a. B ~ R color distribution 
which is rectified for the inclination of the CMR using 
the fit parameters. The galaxies are binned into 0.1 mag 



intervals in the rectified B~R color. An absolute magni- 
tude limit of —20.2 -I- 51ogft.5o (2.0 mag below the mean 
M* in the LF of the cluster galaxies; see Paper III) is 
used to form the distributions. 

The error in the color bin counts is assumed to be Pois- 
sonian. For each cluster a background galaxy color dis- 
tribution is formed by applying the identical rectification 
and binning to the control fields. The background counts 
in color bins are obtained by averaging the control fields, 
and the uncertainty is generated by calculating the stan- 
dard deviation of the counts from the six control fields 
in each bin. The background is then subtracted bin by 
bin from the cluster color distribution, and the total er- 
ror per bin is the sum in quadrature of the error in the 
background and the error in the cluster color bin counts. 

The histogram of the rectified color distribution is 
shown as an insert on the upper left of each CMD in 
Figure 1. On the CMD the -20.2 -|- 5 log /150 mag count- 
ing limit is marked by a dotted line, while the 100% 
completeness magnitude is shown as a solid line. The his- 
togram shows the total count distribution, with the white 
outline representing the contribution of the background 
counts, and the background-subtracted net counts in 
dark shading. 

A Gaussian fit is applied to the resultant color dis- 
tribution to obtain a dispersion. This fit is shown as a 
solid line on the color histogram. Since the color distri- 
bution also contains galaxies other than early-types (as 
in the case when membership and morphological type 
information is available), blue cluster galaxies and clus- 
tered foreground and background galaxies contaminate 
the tails of the distributions. Hence, the Gaussian fit is 
performed by forcing it to fit the peak region properly. 
Visual inspection of the fits indicates that this method 
produces acceptable estimates of the observed width of 
the CMR. 

It is found that the average dispersion of the CMR is 
0.074 mag. The rms of the distribution of the disper- 
sion is 0.026 magnitudes. We can consider this to be 
an upper bound of the true average dispersion of the 
CMR, since a significant amount of the dispersion must 
in fact arise from errors introduced both by photometric 
measurement and stochastic background galaxy contam- 
ination. This is supported by the fact that for clusters 
at lower redshift, where both effects are minimized, the 
CMRs have a lower average dispersion. The mean CMR 
dispersion for the 9 clusters with z < 0.04 is 0.061±0.029, 
where the uncertainty is the rms of the distribution; for 
the remaining 48 clusters, the mean is 0.076 ± 0.025. 
Thus, the mean of the low-redshift sub-sample is smaller 
than the high-redshift sample. 

For comparison, a fainter magnitude cutoff of Mr = 

— 17.5 + 5 log /i5o is used in fitting a Gaussian to the rec- 
tified background-subtracted cluster galaxy counts. For 
a few higher-z clusters whose 100% completeness limit 
(1.0 mag brighter than the 5(t limit) is brighter than 

— 17.5 + 51ogft,5o, we use the completeness limit as the 
counting limit. The mean CMR dispersion obtained us- 
ing this fainter magnitude limit is found to be 0.106 mag, 
with an rms scatter of 0.040 mag. The larger measured 
dispersion of the CMR using a fainter magnitude cutoff is 
a reflection of the increase in the photometric uncertainty 
rather than an increase in the intrinsic CMR dispersion. 
This is supported by the fact that the mean B — R photo- 
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metric error using the —17.5 + 51ogft,5o magnitude limit 
is 0.134±0.164 mag while for the bright magnitude limit, 
—20.2+5 log/i5o, it is 0.040±0.055 mag (the uncertainties 
are the rms of the distribution). 

Estimating the contribution of photometric and back- 
ground contamination errors to the CMR dispersion is 
a complex matter. This is because the error bars are 
not constant over the magnitude range from which we 
compute the CMR dispersion. In a future paper, we will 
analyze the magnitude of the intrinsic dispersion of the 
CMR in detail using Monte Carlo simulations of the data 
set based on measured photometric error distributions. 
For the current paper, we will use a simple estimate to 
show that the correlation between the dispersion width 
and redshift is not due to an evolutionary effect, but re- 
flects the higher photometric uncertainties of the data 
at higher redshifts. Since our color errors are primarily 
dominated by the errors in the B photometry, we can use 
the mean B photometry errors of the galaxies brighter 
than the —20.2 + 51og/i5o mag limit as a proxy for the 
contribution of photometric error to the CMR dispersion. 

In Figure 2 we plot the ratio of the CMR dispersion to 
the average B photometric error as a function of redshift 
of the clusters. From the plot it is clear that at z > 0.04, 
the mean photometric error is ~ 1/3 of the dispersion 
of the CMR; whereas at lower z, the photometric error 
is negligible compared to the CMR dispersion. There 
are a number of clusters showing a significantly larger 
CMR dispersion than the mean photometric error. There 
are 10 clusters with z < 0.04 which have the dispersion- 
to-photometric error ratio larger than ^ 30. These are 
A260, A496, A634, A779, A999, A1142, A2152, A2247, 
A2634, and Coma. If we use the quadratic differences 
of the observed dispersion and the average photometric 
error as a crude estimate to the intrinsic dispersion, we 
find that they have values ranging from 0.028 mag (for 
A634) to as high as 0.126 mag (A2152). We note that 
the cluster with the largest estimated intrinsic disper- 
sion, A2152, belongs to the group of clusters which ap- 
pear to have a contaminating CMR from a higher redshift 
background cluster (see Section 7.2). Hence, it appears 
that both photometric errors and possible contamina- 
tions from background clusters contribute a significant 
part, if not most, of the apparently large dispersion val- 
ues, and we do not find definitive evidence that there is 
a large range of dispersion in the CMRs. 

5.1. Radial Dependence of the CMR 

The measurement of the CMR is based on a region 
surrounding each cluster of 0.2 r2oo in radius. In order to 
look for radial changes in the measured properties of the 
CMR, we have repeated our cluster CMR fits using an 
annulus extending from 0.2-0.5 r2oo- Due to the increase 
in the extent of the measured region, only 33 clusters 
have the necessary areal coverage. A comparison of the 
same 33 clusters for the two different annular regions 
yields a mean CMR dispersion of 0.081 ± 0.030 mag for 
the inner region and 0.090 ± 0.034 mag for the outer ra- 
dial region. The uncertainties are calculated from the 
rms of the distribution. For the slope of the CMR, the 
0-0.2 r2oo region has a mean value of —0.052 ± 0.008, 
while the 0.2-0.5 r2oo outer region has a mean slope of 
—0.053 ± 0.009 (rms uncertainties). Color changes be- 
tween these two regions can be quantified by comparing 



the difference va B — R. Using the CMR fits, the average 
color difference at i? = 17 between the inner and outer 
radial region is A(i3 — i?) = 0.040 ± 0.055, where the un- 
certainty is the rms. The comparison of the mean CMR 
dispersion, slope, and color difference between the inner 
and outer annuli provides evidence of a slight increase in 
the CMR dispersion and a blueward B — R color change 
with increasing cluster-centric radius. This is probably 
due to the increasin g blue fraction with radius (density 
morphology relation; Dressier 1980) and the greater rel- 
ative number of background galaxies (see Paper IV) . 

6. DISCUSSION 

6.1. Elliptical Galaxy Formation and Evolution in 
Clusters 

The classical understanding of the CMR in terms of 
metallicity rests on the idea that massive stellar systems 
are able to experience longer episodes of efficient star for- 
mation. In these systems a large number of supernova 
explosions are necessary to inject enough energy into the 
ISM in order to surpass the binding energy of the galaxy 
to expel the remaining gas a nd halt star form ation activ- 
ity. It has been shown fe.g.. lCarlbergl '1984') that the re- 
moval of enriched gas by supernova-driven winds is more 
efficient in low-mass systems. Hence, giant ellipticals are 
redder due to a longer episode of subsequent star for- 
mation that result in the higher metallicity of their stel- 
lar population; while dwarf ellipticals have had shorter 
episodes of star formation that result in a lower metallic- 
ity of their star population, and hence their bluer colors. 
Passive evolution sets in once the star formation activ- 
ity has ceased. Any recent episode of star formation will 
cause variations of the slope and broadening of the dis- 
persion of the CMR. A possible extreme situation that 
one can consider is that the CMR can break up if in a 
given cluster only the giants or the dwarfs have experi- 
enced some recent episode of star formation. 

At the present epoch the color effects due to dissipa- 
tive galaxy formation cannot be distinguished from age 
effects. However, KA97 have shown that, if the redder 
color of more luminous early- type galaxies are due to age, 
then the CMRs should not exist at high redshift. Ob- 
servational evidence has shown that clusters of galaxies 
with z up to 1.2 have prominent CMRs (see §1). There- 
fore, KA97 concluded that the origin of the CMR is due 
to metallicity effects induced by supernova-driven winds. 
The interpretation that early-type galaxies in clusters 
were formed at redshifts > 2 should produce very simi- 
lar CMRs in current epoch clusters. In the next section, 
we will look at the universality of the properties of the 
CMR and the resultant constraints on the formation and 
evolution of early- type galaxies in clusters. 

6.2. The Universality of the CMR 

6.2.1. The Color of the CMR 

We can compare the rest-band colors of the CMR in 
our sample by computing the k-corrected B — R color of 
the fitted CMR at a fixed absolute magnitude. We pick 
a fiducial absolute magnitude of Mji = —22.2 -I- 5 log /iso , 
the mean M* in the LF of the cluster galaxies (Bark- 
house, Yee, & Lopez-Cruz 2004b, in preparation). The 
B — R color of the red sequence of galaxies in the CMR 
is computed at the apparent magnitude equivalent to 
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the fiducial absolute magnitude. This color is then cor- 
rected to the rest bands by u sing computed k-corrections 
for B and R for E/SO's from lCole man. Wu. fc WeedmanI 
(p980). 

For the 57 clusters with a CMR fit, a mean k-corrected 
B - R color at Mr = M* = -22.2 + 51ogft.5o of 
1.539±0.063 is obtained, where the uncertainty is the rms 
dispersion of the distribution. A conservative estimate of 
the systematic uncertainties in the calibration of the R 
and B photometry is ^ 0.03 mag for each filter, pro- 
ducing a systematic uncertainty of 0.042 mag m. B — R. 
Hence, the intrinsic dispersion of the color of the CMR 
in this large sample over a small redshift range is about 
0.05 mag. This small dispersion is very similar to that 
estimated by Small et al. (1998), who, using 10 clusters 
at 0.22 < z < 0.28, found also a small observed disper- 
sion of the color of the CMR of ~ 0.04 mag, consistent 
with an intrinsic dispersion of 0.03 mag. Our measured 
int rinsic scatter is also consistent with value s measured 
bv lStanford. Eisenhardt. fc Dickinsm] l|1998j) for a sam- 
ple of 19 clusters spanning a redshift range from 0.03 to 
0.9 (see their Fig. 6). 

6.2.2. The Dispersion of the CMR 

The dispersion of the CMR can potentially serve as 
a statistical indicator of the star formation history of 
early- type cluster galaxies (e.g., BLE). For example, if 
the dispersion is connected to the initial formation, it 
could be used to estimate the spread of the formation age 
of the galaxies. Or alternatively, the dispersion could be 
an indication of additional episodes of star formation. 

Typically, most investigations have found relatively 
small dispersions in the CMRs, although they are all 
based on very small samples and usually sampling only 
over a range of four or five magnitudes. For example, 
BLE found a very small dispersion of 0.05 mag inU — V 
with about 0.03 attributable to photometric errors. In 
B — R, and over a similar baseline of magnitude, we 
obtain a comparable measurement for Coma — a dis- 
persion of 0.06 mag with a mean photometric error of 
0.01 mag, using the background subtraction technique 
(with no membership and morphological type informa- 
tion). Using HST data with redshift a nd mo rphological 
type information, Ivan Dokkum et al] l)1998fl performed 
very detailed analysis of the CMR of the rich EMSS clus- 
ter 1358-1-62 at z = 0.33, and found again very small 
dispersions of 0.02 to 0.04 mag, depending on the ex- 
act galaxy sample chosen. .Ellis et al,. (,1997.) , using HST 
data for the core regions of three z ~ 0.55 clusters, found 
dispersions ranging from 0.07 to 0.11 mag. F o r high red- 
shift clusters {z^ 1. 3), Ivan Dokkum et aP l)200lD and 
iBlakeslee et"al] (|2^, using HST data, measured dis- 
persions of 0.04 and 0.03 mag, respectively. Our mean 
dispersion of 0.074 mag can only serve as an upper limit 
to the ensemble average, and is consistent with these 
measurements. Despite our large sample, we have found 
no definitive evidence of clusters with dispersion broader 
than those reported in the literature, i.e., larger than 
~ 0.1 mag. This lack of broad dispersion clusters is con- 
sistent with an early formation epoch and the lack of 
significant recent star formation in the early-type galax- 
ies in clusters. 

6.2.3. The Slope of the CMR 



Comparing the variation of the slope of the CMR to 
those predicted by models is more robust than com- 
parisons based on absolute color evolution. Uncertain- 
ties in the transformations to standard photometric sys- 
tems, both in the observations and the models, and k- 
corrections can considerably reduce the significance of 
any photometric test; whereas t he slope of the CMR is 
subject only to internal errors. IGladders et all llT99l 
exploit this robustness of the comparison between mod- 
els and observations of the slope of the CMR for clus- 
ters at z ~ to 0.8. For the low-redshift anchor, they 
used redshift-binned composites from a sub-sample of 
the full clus ter sample CMRs that were available from 
lL6Dez-Cru25 (1997). Comparing with CMRs of higher 
redshift clusters from HST data, IGladders et alJ fl998') 
concluded that the formation epoch, assuming a coeval 
early-type population in clusters, is conservatively larger 
than z ~ 2.5. 

The 57 clusters from our sample provides one of the 
largest single samples of CMRs over a small redshift 
range obtained in a homogeneous fashion. This allows 
us to compare a large sample of clusters of various rich- 
ness and cluster morphological types, at about the same 
epoch, to test the universality of the CMR in a single 
epoch. In Figure 3 we plot the distribution of the B — R 
CMR slopes of individual clusters as a function of the 
cluster redshift. The slopes have a very narrow range of 
values between —0.03 to —0.08, with a definite trend of 
steeper values towards larger redshifts. 

To interpret the slopes of the CMRs, we appeal to the 
models of KA97. The KA97 models are sophisticated 
stellar population synthesis models which include new 
libraries for stellar evolution accounting for the effects of 
metallicity. Kodama (1996, private communication) has 
repeated the calculations in KA97 for the filters employed 
by this study. The transmission curves for the B, R 
and / filters were obtained from KPNO. An epoch of 
galaxy formation of 15 Gyr (zf « 5.4) is assumed. A 
more critical parameter is the epoch when the supernova 
driven winds {tguj) are established. KA97 have calibrated 
igio using the observed CMR for Coma (BLE) based on 
U — V colors. The resultant variations of the model slope 
as a function of redshift is plotted on Figure 3 as a dashed 
line. We remark that this model is not a fit, i.e., no 
freely adjustable parameter based on our cluster data has 
been applied. The model was generated independently 
using as the anchoring calibration BLE's data for Coma. 
The model clearly provides an excellent fit, including the 
slight steepening trend of the slope with redshift. The 
very slight apparent slope variations seen here are due 
to the CMRs being sampled at slightly bluer rest wave- 
bands at higher redshifts, and it is not due to stellar 
evolutionary effects. Hence, once the CMR is properly 
corrected for this band pass effect, we can conclude that 
the CMR is very uniform in its parametrization for the 
whole sample of 57 clusters. 

6.2.4. The Universality of the CMR and Cluster 
Evolution 

Our cluster sample contains a diverse set of clus- 
ters: it covers a richn ess range over a factor of 4 
ijYee fc L6pez-Cni3ll999 ^. from Abell to amongst the 
richest known clusters; it contains clusters that cover the 
full range of the morphological types for schemes in the 
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optical and in the X-rays fsee lLopez-CrudEoOSl for a re- 
view on the morphological classification of clusters). The 
lack of variations in the slope of the CMR, along with the 
very small dispersion in the color, across all clusters with 
such diverse properties allows us to draw the very conser- 
vative conclusion that no cluster has experienced a sig- 
nificant star formation episode in the core during the last 
3 Gyr (i.e., the look-back time of z ~ 0.2). Furthermore, 
this universality indicates that all current epoch clus- 
ters contain a core set of early-type galaxies which were 
fo rmed at a suffic i ently distant epoch. As demonstrated 
in lGladders et alJl)1998t) . using models from KA97, it can 
be shown that a significant digression of the CMR slope, 
in the form of a turn-over to a shallower slope, is expected 
due to the differing evolution rates of early-type galaxies 
with different metallicities when the galaxies are at ages 
younger than ~ 4 Gyr after formation. Hence, we can 
conclude that none of these 57 clusters were formed less 
than 7 Gyr ago (or z ~ 0.8 for Hq = 70 km s~^ Mpc~^, 
rtm — 0.3, and SIa — 0.7). Given the relatively large 
sample size, one can conclude that it is very unlikely 
that there are any current-epoch clusters (rich enough 
to be in the Abell sample) which w ere formed le ss than 
7 Gyr ago. We note that fdaddcrs et aP l|1998() . com- 
bining 44 clusters from the LOCOS sample and HST 
data for 6 higher redshift clusters, suggested, based on 
the variation of the CMR slope with redshift, that early- 
type galaxies in clusters form an old population with for- 
mation redshift Zf > 3. However, because of the small 
number of high-z clusters and their high average rich- 
ness, they could not rule out that some (likely less mas- 
siv e) clusters may form at low er re dshifts. Recent s tudies 
bv lvan Dokkum et alJ l|200lD and lBlakeslee et alJ l)2003l^ 
extend the redshift limit to z ~ 1.3 for the existence of 
a well-established early-type cluster population, further 
supporting a formation redshift z/ > 3. A large sample 
of clusters covering a wide range of cluster properties at 
z between 0.5 and 1.0 will allow us to place stringent con- 
straints on the formation epoch of clusters as a function 
of cluster properties. 

6.3. The Coma Cluster and the CMR for Dwarf 
Galaxies 

The Coma cluster (Abell 1656) is the most studied 
galaxy cluster; its proximity and richness makes it ideal 
for all varieties of cluster studies. It was included as an 
anchoring reference for the properties that we measure 
in the other clusters of our sample. What is striking is 
that, because of the low redshift of the cluster and the 
depth of sampling, its CMR (see Figure 1) can be defined 
over a range of 8 magnitudes {R — 12.5 to 20.5) easily, 
and it can be detected well into our 100% completeness 
limit of i? = 21.2 (~ —14.2 -|- 51og/i5o in absolute mag- 
nitude). If we de fine the dwarf galaxy popul a tion i n the 
same manner as iSecker. Harris, fc Plummeij l)1998f) . i.e, 
galaxies fainter than i? = 18 in apparent magnitude (ab- 
solute magnitudes fainter than —18-1- 51og/i5o), it is ap- 
parent that the same linear fit describes the CM R from 
giants to dwarfs. This result agrees withj^andagc (|1972fl 
indication that dwarf galaxies follow the CMR defined 
by bright elliptical galaxies in Virgo. Visual inspection 
of the objects on the sequence zone indicates that their 
morphology resembles dwarf elliptical galaxies. Hence, 
we do not detect any change in slope in the CMR over 



this very large range of galaxy masses. We also note that 
we do not detect any significant changes of slope within 
the CMR in any of our clusters, most of which we can 
sample to and past absolute magnitude — 17-f 51og/i5o 
(although, because of the higher redshift, the background 
confusion is usually larger). Thus, if there is a bend in 
the CMR, as has been claimed by a number of studies 
(e.g., Metcalfe et al. 1994), this is not a frequent oc- 
currence at these redshifts. The bending may be due to 
photometric errors or due to the working definition of 
the color aperture. 

The slope of the CMR can be used a s an age in- 
dicator (§6.2.3 and IGladders et all I1998|) : hence, the 
fact that dwarfs appear to follow the same CMR as 
the giants suggests that dwarf galaxies have a simi- 
lar evolutionary history as giant ellipticals. This seem 
to be the case: lOdell. Schombert. fc Rako'^ l)2002|) us- 
ing medium-band photometry in the Coma cluster have 
found that dwarf galaxies do follow the CMR defined 
by gia nt carly-typc galaxies. However, according to 
lOdell. S chombert. & Rakps dwarf ellipticals, depending 
on the morphology, could be either younger or older than 
giant ellipticals. 

We have explored some observational properties of the 
CMR. We found that these properties can be explained 
within the classical models. However, alternative mod- 
els can be intro duced considering both age and metallic- 
ity effects fe.g.. lFerreras. Chariot^ fc Silklll999D or alter- 
nativ e formation scenarios (e.e.. IKauffmann fc Charloti 
119981) . A complete scenario for the formation of early- 
type galaxies in clusters should also incorporate other 
observational constraints, i.e., line strength indices and 
their evolution, and the fundamental plane as its evolu- 
tion. 

7. PRACTICAL APPLICATION OF THE CMR 

The CMR, beside being an important diagnostic in 
the formation and evolution of early- type galaxies, can 
also serve as a powerful observational tool for detecting 
galaxy clusters and measuring photometric redshifts. In 
the following we briefly examine these applications. 

7.1. A Single- Color Redshift Estimator 

In the last few years, there has been a resurgence of in- 
terest in and applications of photometric techniques for 
determining redshifts of galaxies (sec, e.g., Ycc 1999, for a 
review). Typically such techniques use about four filters, 
combined with either empirical training-set calibrations 
or model spectral energy distributions, to derive the red- 
shift. The CMR of clusters of galaxies, however, offers 
the unique capability of serving as a redshift estimator 
using only two filters, since the population type of the 
galaxies is known. In the observer's frame the CMR of 
a cluster shifts toward redder colors as one progresses in 
redshift. It is possible to calibrate the changes in color of 
the CMR using our cluster sample as a training set and 
generate a redshift estimator. This indicator, calibrated 
empirically in the observer's frame, will also account for 
evolutionary effects, if any is detectable at these low z. 
In the following we determine an empirical correlation 
between the B — R color of the CMR and redshift. 

The zero points from the fits of the CMR (which is 
equivalent to the color of a galaxy at apparent magni- 
tude) in Table 1 show a strong correlation with redshift, 
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but the errors are large, as this value is greatly corre- 
lated with the uncertainty in the slope. A much bet- 
ter estimator can be obtained by using a point on the 
CMR within the magnitude range of the observed CMR. 
We have chosen a fiducial, but arbitrary, magnitude of 
R = 17. An apparent magnitude is chosen as the fidu- 
cial because no prior knowledge of the cluster's redshift 
is assumed when determining the redshift of an unknown 
cluster. The {B — R)ji^n values are listed in column 6 of 
Table 1. The average error in {B — R)b.=i7 as estimated 
from the fit to the CMR by bootstrap is ~ 0.016 mag 
for clusters at z 0.04, and 0.030 mag for clusters at 
z ~ 0.15. In Figure 4 we plot the very tight correlation 
between {B — R)u=i7 with redshift. The correlation be- 
tween the computed (B — R)fi^n and z is fitted with a 
quadratic polynomial anchored on Coma: 

z = 0.0849(+/ - 0.0165)(B - Rf - 0.1419(+/ - 0.0573) 
(B-i?) + 0.0511(+/- 0.0496). (2) 
This fit is shown as a solid line in Figure 4. The disper- 
sion about the fit is cr = 0.010, or a mean Az/z ^ 0.13. 
This dispersion is comparable or better than the best 
az ^ 0.04 that is typically obtained using training set 
techniques for redshift determination using color infor- 
mation from at least fou r bands for z ~ 0.3 galaxies 
fe.g.. 'Brunncr et alJll99'^ . The high accuracy of this 
method rest in the fact that we are using a large number 
of galaxies with a known spectral type simultaneously. 
iPostman fc Laueil l) 19951) have devised a redshift indica- 
tor based on the Lm — cx (the metric luminosity and the 
logarithmic slope of the surface brightness profile) rela- 
tionship for BCGs. For our sample, we found that the 
mean and standard deviations for (^°'"~^"') were -0.0164 
and 0.180, respectivel y. These values are comparable to 
the ones found in the 'Post man fc LaueJ study; however 
we have covered a redshift range about four times as 
large. It is interesting to note that the CMD of A2152 
(Figure 1) shows the presence of a contaminating clus- 
ter at B — R ^ 2 , or z = 0.11 according to equation 
1. iBlakesle^ ()200lD has recently reported the discovery 
of a cluster 2.4' from the center of A2152 at a spectro- 
scopically determined redshift of z — 0.13, and having 
a red sequence at B — R ^ 2.1. The CMR provides a 
very reliable redshift estimator that is based on a uni- 
form property of elliptical galaxies. 

7.2. The CMR as a Tool for Cluster Detection 

The universality of the CMR, at least at the low- 
redshift range being considered, suggests that the sig- 
nature of a CMR in a CMD is a good indicator of galaxy 
clusters. In the case of our data, such a signature can 
be used to detect high-density concentrations of galax- 
ies in the background or foreground of the target Abell 
clusters. Simple visual inspection of the CMDs (Figure 
1, and Paper I) shows that the presence of such contam- 
inating clusters is very easily seen. Prominent examples 
include A690, A2152, and A2399. However, a more ef- 
ficient way to detect a background cluster is to use the 
rectified color histograms that are applied to determine 
the dispersion of the CMR (Section 5). A background 
cluster can be easily seen, after background subtraction, 
as a strong peak with a median color redder than the 
foreground target cluster. By inspection of the color his- 
tograms, we identify nine clusters (A514, A690, A999, 



A1291, A2152, A2247, A2271, A2399, and A2634) as 
having a very high probability of having a background 
cluster within their image frames, and an additional four 
(A168, A671, A1795, and A2440) as having a possible 
background cluster superposition. A robust and rigor- 
ous cluster-finding algorithm based on this idea of the 
CMR as a cluster detection signature, which searches for 
enhanced density in the four-dimensional space of color, 
magnitude, and x-y sk y position has been developed by 
'Glad ders fc Yed 1)2000(1 . and this algorithm will be ap- 
plied to our cluster data to quantify the background clus- 
ter detections in a future paper. 

At the modest depth of our survey, which allows us to 
detect clusters out to about z ^ 0.3, using this crude 
method of identifying clusters from the B — R color dis- 
tribution, we have detected 9 to 13 background clusters 
to our target sample of 57 clusters, or at a rate of about 
15 to 20%. The B — R colors of the background clusters 
are typically greater than 2.3 mag, indicating they are at 
a redshift of greater than 0.15 (see Section 7.1). In total, 
we have 57 cluster observations sampling ^ 15 square 
degrees, with an average completeness limit of i? = 21.2. 
Thus we detect ~ 1 cluster per square degree having a 
median redshift of z = 0.2. 

There is a general perception that X-ray surveys and 
observations are free from background cluster-cluster 
contaminations. It is interesting to examine the few 
cases where we find definite optical projections. A 
remarkable case is A514, which has been considered 
standard cluster for substructure studies in the 
X-r ay (e.g., Ebcling ct al. 1996; J ones & Forman 199i 
iGovoni et alJl200lt iKolokotronis et all 120011) . The X- 
ray image of this cluster shows many substructures, 
including a very pro minent component to the NW 
l| Jones fc Formanll999j) . The CMD clearly shows a back- 
ground CMR at the color of B — i? ~ 2.8, correspond- 
ing to a redshift of greater than 0.20, and the galax- 
ies contributing to this background CMR are coincident 
with this large substructure. Thus, at least part of the 
complex X-ray structure of A514 arises from the super- 
position of clusters. Another very clear case of cluster 
superposition in the X-ray is provided by A2152. The 
X-ra y image shows a strong source with a NW exten- 
sion l)Jones fc Form an 1999). In this case, the CMD 
clearly shows the presence of a background cluster at 
i? — i? ~ 2.0 and has been spect roscopically verified as 
a z = 0.13 cluster (Blakeslee 2001), which is coincident 
with the X-ray extension. Of the remaining seven cases 
of clear background CMRs, only A999 and A1291 do not 
show significant X-ray substructure. Furthermore, of the 
four clusters with possible superposing CMDs, two show 
significant X-ray substructure. Hence, we have found 
that in a number of instances, background cluster-cluster 
superpositions may well have a significant effect on the 
apparent X-ray morphology of clusters. A more detailed 
analysis of the cluster-cluster superpositions to study the 
correlation between such superpositions and X-ray struc- 
tures based on the more quantitative red - sequen ce clus- 
ter finding algorithm of iGladders fc Ye3 l|200(]f) will be 
presented in a future paper. 

8. CONCLUSION 

In this paper we have explored the properties and ap- 
plications of the CMR using a large sample of 57 low- 
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redshift Abell clusters. We demonstrate that simple tech- 
niques based on the biweight function are able to fit a sin- 
gle line to the CMR. A straight line describes the CMR 
for a range of at least 5 magnitudes for all the clusters in 
the sample, and extending as far as 8 magnitudes when 
the data are of sufficient depth. We find no evidence for 
a change of slope within the CMR over this magnitude 
range. 

The properties of the CMR for this diverse sample 
of low-redshift clusters are remarkably universal. We 
also find no definitive evidence that the dispersion of 
the CMR is significantly greater than about 0.1 mag 
in any cluster. The distribution of the k-corrected 
color at a fixed absolute magnitude of the individual 
CMRs in the sample have a very small dispersion of 
0.05 mag. Finally, the slope of the CMR is entirely 
consistent with being identical, once the CMRs are cor- 
rected to the same rest wave bands. Given the con- 
stancy of the parameters of the CMR, the universality 
of this fundamental property of early-type galaxies in 
clusters is established. In agreement with the results of 
IStanford. Eisenhardt. fc DickinsoiJ l)1998f) . we show that 
our observations are consistent with models that explain 
the origin of the CMR based on galaxy formation which 
produces the color dependence due to varying metal- 
licities. The universality of the slope and color of the 
CMR enable us to place a conservative limit that the last 
episode of significant star formation in cluster early-type 
galaxies happened more than ^ 3 Gyr ago, and that the 
core set of early-type galaxies in current epoch clusters 
were formed more than ~ 7 Gyr ago. 

We demonstrate that the CMR provides a powerful 
tool for the detection of clusters and also can be used as 
a very accurate and efficient redshift indicator for groups 
and clusters of galaxies. The CMR can also be effec- 
tively used to limit the contribution of background galax- 
ies in the statistical determination of cluster member- 
ship, in that galaxies redder than the CMR can be con- 
fiden t ly considered as background galaxies llLopez-Crud 
ri997tlL6nez-Cniz et abllTflflTt ma,rkhousell2nn3D. We will 
exploit this property in the companion paper on the lu- 
minosity function of cluster galaxies (Paper III). 
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Fig. 1. — Example color-magnitude diagrams from our cluster sample in B — ii vs R for the total region covered by our images. The solid 
line shows the 100% completeness limit, while the dotted line indicates absolute magnitude —20.2 at the cluster's redshift. The dashed line 
illustrates the best fitting color-magnitude relation for the cluster, with the fit equation shown at the lower left. The insert histogram is 
the OMR rectified B — R color distribution for background-subtracted galaxies within 0.2 r2oo of the adopted cluster center. The Gaussian 
fit to the OMR is shown as the thick solid line over-plotting the histogram with the dispersion indicated on the upper left corner. A total 
of 6 CMDs are shown. These are: (a) A21, (b) A690, (c) A1213, (d) A1656 (Coma), (e) A2152, and (f) A2399. The CMD of the complete 
sample can be found in Paper I. We note that A690, A2152, and A2399 are examples of cluster fields which have an apparent ba<;kground 
CMR, likely arising from background clusters at z > 0.2. These can be detected both in the CMDs and as the prominent second peak in 
the B — R histograms. 

Fig. 2. — The ratio of the computed dispersion of the CMR from the B — R histogram to that of the mean B-band photometric uncertainty 
of the galaxies versus redshift. Note that for clusters with z > 0.04, their CMR dispersion is partly due to photometric errors. For the 
z < 0.04 clusters, there arc a number of clusters which show a significantly larger CMR dispersion than expected simply from measurement 
errors. At least part of this dispersion can be attributed to excess background contaminations from cluster-cluster superpositions (see text). 

Fig. 3. — The variation of the slope of the CMR with redshift. The points correspond to the cluster sample. The dashed line is not a fit, 
but the predicted variation derived from a spectral synthesis model that accounts for variations due the metallicity and passive evolution 
in elliptical galaxies (K A97) . The equation parameterizing this model as a function of z is given at the top of the graph. 

Fig. 4. — The variation of the CMR color at a fixed apparent magnitude with redshift. The points correspond to the B — R colors at 
R = 17.0 for the clusters in this survey as computed from the best fitting CMR. The solid line represents a second order polynomial fit to 
the data (see text for details). 
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TABLE 1 

Cluster Sample and CMR Fits 



Cluster 


Redshift(2) 


Slope (b) 


Zero Point (a) 


^CMR 


(B - R)R: 


A21 


0.0946 


-0.057 


-I- 


0.012 


2.88 


-I- 


0.20 


0.080 


1.91 


A84 


0.1030 


-0.052 


± 


0.013 


2.75 




0.21 


0.097 


1.87 


A85 


0.0556 


-0.038 


-I- 


0.014 


2.31 




0.22 


0.113 


1.66 


A98 


0.1043 


-0.056 




0.013 


2.99 




0.22 


0.077 


2.04 


A154 


0.0638 


-0.059 


-I- 


0.013 


2.77 


± 


0.22 


0.115 


1.77 


A168 


0.0452 


-0.048 


-j- 


0.013 


2.34 


± 


0.20 


0.042 


1.53 


A260 


0.0363 


-0.045 


± 


0.014 


2.30 


± 


0.22 


0.093 


1.53 


A399 


0.0715 


-0.048 


-I- 


0.013 


2.58 


± 


0.22 


0.088 


1.76 


A401 


0.0748 


-0.047 


-I- 


0.012 


2.66 


± 


0.19 


0.062 


1.86 


A407 


0.0463 


-0.046 


-I- 


0.013 


2.47 


± 


0.21 


0.072 


1.69 


A415 


0.0788 


-0.044 


-I- 


0.013 


2.74 


± 


0.23 


0.055 


1.98 


A496 


0.0329 


-0.052 


-I- 


0.013 


2.43 


± 


0.22 


0.063 


1.54 


A514 


0.0731 


-0.048 


-|- 


0.013 


2.56 


-£ 


0.21 


0.075 


1.75 


A629 


0.1380 


-0.050 


-I- 


0.010 


3.00 


-I- 


0.19 


0.112 


2.15 


A634 


0.0265 


-0.047 


-I- 


0.014 


2.40 


-1- 


0.22 


0.028 


1.61 


A646 


0.1303 


-0.058 


-|- 


0.010 


3.10 


-|- 


0.18 


0.130 


2.11 


A665 


0.1819 


-0.076 


-j- 


0.012 


3.56 


-|- 


0.21 


0.145 


2.27 


A671 


0.0502 


-0.045 


-I- 


0.014 


2.41 


-I- 


0.22 


0.070 


1.64 


A690 


0.0788 


-0.041 


-)- 


0.011 


2.52 


-£ 


0.19 


0.091 


1.83 


A779 


0.0229 


-0.047 


-(- 


0.016 


2.34 




0.24 


0.109 


1.54 


A957 


0.0450 


-0.049 


-I- 


0.013 


2.44 


± 


0.21 


0.061 


1.60 


A999 


0.0323 


-0.054 


± 


0.017 


2.46 


± 


0.24 


0.085 


1.54 


A1142 


0.0349 


-0.031 




0.013 


2.09 


± 


0.22 


0.034 


1.56 


A1213 


0.0468 


-0.043 


-j- 


0.012 


2.34 




0.19 


0.073 


1.60 


A1291 


0.0530 


-0.050 


-I- 


0.012 


2.58 


± 


0.21 


0.077 


1.72 


A1413 


0.1427 


-0.056 


-I- 


0.013 


3.08 




0.23 


0.076 


2.13 


A1569 


0.0784 


-0.064 


± 


0.013 


2.86 




0.22 


0.067 


1.77 


A1650 


0.0845 


-0.052 


-I- 


0.012 


2.70 


± 


0.22 


0.084 


1.82 


A1656 


0.0232 


-0.046 


-I- 


0.012 


2.22 




0.21 


0.056 


1.44 


A1775 


0.0717 


-0.050 


-I- 


0.014 


2.58 




0.24 


0.071 


1.74 


A1795 


0.0622 


-0.059 


± 


0.013 


2.77 


± 


0.23 


0.067 


1.77 


A1913 


0.0528 


-0.050 


± 


0.013 


2.44 


-I- 


0.21 


0.047 


1.60 


A1983 


0.0449 


-0.048 




0.013 


2.40 


-I- 


0.22 


0.062 


1.59 


A2022 


0.0578 


-0.048 


-£ 


0.013 


2.45 


-1- 


0.22 


0.070 


1.64 


A2029 


0.0768 


-0.050 


± 


0.012 


2.67 


-I- 


0.20 


0.080 


1.82 


A2152 


0.0410 


-0.044 


± 


0.011 


2.38 


-I- 


0.20 


0.126 


1.63 


A2244 


0.0968 


-0.061 


-J- 


0.012 


2.96 


-|- 


0.21 


0.073 


1.92 


A2247 


0.0385 


-0.060 


-)- 


0.011 


2.50 


-(- 


0.18 


0.032 


1.48 


A2255 


0.0808 


-0.046 


± 


0.012 


2.59 


± 


0.20 


0.073 


1.81 


A2256 


0.0581 


-0.047 


-|- 


0.013 


2.46 




0.22 


0.050 


1.67 


A2271 


0.0576 


-0.051 


-I- 


0.012 


2.72 


± 


0.20 


0.059 


1.85 


A2328 


0.1466 


-0.063 


-|- 


0.014 


3.14 


-£ 


0.23 


0.130 


2.06 


A2356 


0.1161 


-0.060 


-|- 


0.013 


3.07 


-£ 


0.24 


0.061 


2.05 


A2384 


0.0943 


-0.059 


-1- 


0.013 


2.93 


-£ 


0.23 


0.046 


1.93 


A2399 


0.0587 


-0.046 


± 


0.013 


2.52 


± 


0.22 


0.070 


1.74 


A2410 


0.0806 


-0.044 


± 


0.014 


2.65 


± 


0.24 


0.031 


1.90 


A2415 


0.0597 


-0.044 


± 


0.013 


2.48 


± 


0.23 


0.089 


1.73 


A2420 


0.0838 


-0.048 


± 


0.013 


2.75 


± 


0.23 


0.078 


1.93 


A2440 


0.0904 


-0.057 


± 


0.013 


2.96 


± 


0.23 


0.088 


2.00 


A2554 


0.1060 


-0.051 


± 


0.014 


2.92 


± 


0.23 


0.052 


2.05 


A2556 


0.0865 


-0.059 


± 


0.012 


2.94 


± 


0.21 


0.049 


1.93 


A2593 


0.0421 


-0.041 


± 


0.014 


2.34 


± 


0.23 


0.048 


1.65 


A2597 


0.0852 


-0.070 


± 


0.012 


3.12 


± 


0.24 


0.037 


1.94 


A2626 


0.0573 


-0.043 


± 


0.012 


2.51 


± 


0.20 


0.059 


1.78 


A2634 


0.0310 


-0.053 


± 


0.013 


2.40 


± 


0.20 


0.047 


1.50 


A2657 


0.0414 


-0.054 


± 


0.013 


2.68 


± 


0.21 


0.077 


1.75 


A2670 


0.0761 


-0.046 


± 


0.012 


2.64 


± 


0.20 


0.094 


1.86 



